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Abstract. We propose a method for recovering the two-dimensional gravitational potential of galaxy clusters which combines 
data from weak and strong gravitational lensing. A first estimate of the potential from weak lensing is improved at the ap- 
proximate locations of critical curves. The method can be fully linearised and does not rely on the existence and identification 
of multiple images. We use simulations to show that it recovers the surface-mass density profiles and distributions very accu- 
rately,even if critical curves are only partially known and if their location is realistically uncertain. We further describe how 
arcs at different redshifts can be combined, and how deviations from weak lensing can be included. 



1. Introduction 

Weak gravitational lensing constrains the projected mass dis- 
tributions of galaxy clusters with an angular resolution of 
> 0.5', while strong lensing occurs typically not farther than 
(0.5 - 1)' from cluster cores. Nonetheless, both phenom- 
ena are due to the same gravitational potential. Recent ob- 
servations of galaxv clusters with the ACS camera on-board 
HST (see iBroadhurst et alJ l2005l for an example) have re- 
vealed large numbers of arcs in individual galaxy clusters, for 
which Abell 1689 is an outstanding example. In such clus- 
ters, strongly lensed images provide numerous constraints on 
the lensing potential, and the question is raised how weak and 
strong lensing are best combined in joint reconstructions of the 
lensing mass distribution. 

Several methods were recently proposed whi ch rely on 
the identification of multi ply-imaged systems (Bradac et alJ 
l2005albtlDiego et alJ2005l) . We propose an alternative method 
here for which multiple images are not necessary and thus 
do not need to be identified. Based on a least-;^'^ minimi- 
sation, a first estimate of the lensing potential is obtained 
from weak lensing alone, following the approach suggested 
bv Bartelmann et al. (1996) and Se itz et a l. (1998). Where the 
position of critical curves can be approximately identified due 
to the presence of strongly-lensed images, this estimate is im- 
proved by requiring that the determinant of the Jacobian matrix 
vanish. The method is easily generalised to sources at dififerent 
redshifts. It can be fully linearised and allows deviations from 
the strict weak-Iensing limit to be accounted for by simple iter- 
ation. 

We outline the main ideas underlying the method in Sect. 2. 
After reviewing the necessary formalism of gravitational lens- 



ing in Sect. 3, we describe the method in Sect. 4 and il- 
lustrate its feasibility using numerical simulations in Sect. 5. 
Conclusions are given in Sect. 6, and technical detail is sum- 
marised in the Appendix. 



2. Overview 

The method we are proposing rests on two central ideas. First, 
both weak and strong lensing must be described by the same 
underlying, two-dimensional lensing potential. Data from weak 
and strong lensing observations must thus be expressed in 
terms of constraints on the potential. Second, strong lensing 
can be located in several ways, of which the identification of 
multiple images is one, and the existence of critical curves is 
another Recent images of strongly-Iensing galaxy clusters il- 
lustrate that gravitational arcs in many cases allow an accurate, 
piece-wise locaUsation of critical curves. On the other hand, 
multiple images need to be identified first and are only caused 
by sources close to and inside caustics, while sources outside 
caustics can still be used to identify critical curves. 

Based on these ideas, our algorithm can be summarised 
very simply: we wish to find a map of the lensing potential, 
discretised on a grid, which is determined such as to optimise 
the agreement with the observed data. The gravitational shear 
caused by weak lensing imprints coherent distortions on the 
images of background galaxies. To lowest-order approxima- 
tion, their average ellipticity is an unbiased estimator of the 
local shear The lensing potential must be aiTanged such as 
to reproduce this measured shear. Critical curves occur at sin- 
gular points of the lens mapping, which locally constrain the 
curvature matrix of the potential. Thus, our algorithm aims at 
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optimally constraining the potential such that its tidal field re- 
produces the measured shear, and its curvature reproduces the 
critical curves. 

Several complications come to mind immediately. First, 
there is a problem of scales. Given number densities of ~ 
30 - 40 background galaxies per square arc minute, the angu- 
lar resolution of weak-lensing observations is limited to > 0.5'. 
Strong lensing, however, typically happens within < 1' of clus- 
ter centres. Both can be combined only if grids of high reso- 
lution in the core and low resolution farther outside are used. 
For the practical implementation, we choose to first reconstruct 
the potential at low resolution on a coarse grid, then refine the 
grid everywhere in the core area where critical curves exist, fill 
the refined grid by suitable interpolation from the coarse grid, 
and correct the potential there such as to reproduce the critical 
curves. 

Second, galaxy ellipticities reflect the shear y only in the 
limit of weak lensing, a- <k 1 . Deviations occur towards criti- 
cal curves. Outside critical curves, ellipticities constrain the re- 
duced shear y/(l - k), and its inverse comp l ex conj ugate 1 /g* 
inside. As other studies (e.g. Brada c et alJ ll2005b)') have re- 
marked earlier, this is not a problem of principle because it can 
be solved in a quickly converging iteration. 

Third, arcs in clusters typically originate from sources at 
diflerent redshifts. Since the geometrical efficiency of a fixed 
lens increases monotonically with source redshift, so does the 
lensing potential. Arcs at different redshifts thus constrain the 
curvature of a lensing potential with the same shape but differ- 
ent amplitude. However, the potential grows with source red- 
shift by a linear factor This suggests the introduction of a ref- 
erence potential for one arbitrary source redshift, to which the 
potentials for sources at other redshifts can linearly be trans- 
formed. 

Fourth, critical curves can almost never be traced through- 
out a cluster, but only piece-wise. This is irrelevant for our pur- 
poses since the constraint that the curvature of the lensing po- 
tential must reproduce the critical curves is purely local. As the 
distortion of arcs close to critical curves varies rapidly with po- 
sition, measurement errors for the location of critical curves are 
typically small. Alternatively, constraints on critical curves ob- 
tained from parameterised lens models can be used in a hybrid 
approach. 

Finally, all constraints can be expressed in terms of second 
derivatives of the lensing potential. In the limit of weak lens- 
ing, the constraints are linear since the source ellipticity is then 
a linear combination of second potential derivatives. Moving 
into the strong-lensing regime, the ellipticity constraints can 
be kept linear by means of an iterative procedure. We devise a 
scheme in which even the constraints from critical curves can 
be expressed by a linear correction term. While this linearity 
is not conceptually important, it is practically because it allows 
minimisations by matrix inversions. 

3. Basic formalism 

We start by reviewing the basic formalism for gravitational 
lensing as we shall need it in the course of the paper. Dealing 
with isolated lensing systems such as galaxy clusters, we adopt 



the thin-lens approximation, according to which the lensing 
mass distribution is projected onto a lens plane perpendicular 
to the line-of-sight. Sources are located on source planes which 
are also perpendicular to the line-of-sight. The lens system is 
characterised by the three angular-diameter distances Di_s,is be- 
tween the observer and the lens, the observer and the source, 
and between lens and source, respectively. 

All relevant properties of the lens system are then contained 
in the scalar lensing potential which is the suitably projected 
and rescaled Newtonian gravitational potential of the lens. 
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(e.g. ISchneideret all ll 9921) . Obviously, the lensing potential 
depends on the source redshift. Assuming that the lensing mass 
distribution is the same for sources on different source planes, 
we can still introduce a single scalar potential \p{Q) for a fiducial 
source redshift z^, and then scale the potential to other source 
redshifts Zs according to 



Du(zs) DsiW 



(2) 



The linearity of this transformation allows the linear recon- 
struction of the single, fiducial lensing potential ^ from sources 
on multiple source planes. 

The two-dimensional, projected mass distribution of the 
lens is described by the convergence k, which is the surface 
mass density S in units of 
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Distortions caused by the lens are characterised by the trace- 
free, symmetric shear tensor with the two components yi_2. 
Both are linear combinations of second derivatives of i//. 
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7i = :^i4f^n -(1^,22) , 72 -^,12, k - -V ^ . (4) 

The lens mapping relates the source position fi to the image 
position(s) 6, 



(5) 



For sources which are small compared to typical scales 
of the (reduced) deflection angle a = Viff, the lens mapping 
can be linearised. Imaging is then locally characterised by the 
Jacobian matrix 
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This shows that k is responsible for shrinking or stretching im- 
ages isotropically, while y causes anisotropic deformation. 

A sufficiently small circular source of radius r is imaged as 
an ellipse with semi-major and -minor axes a = r(l - k - y)"' 
and b - r(l - K + y) ', respectively. The ellipticity, defined as 
e = (a - b)/{a + b), is thus 



e = 
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defining the so-called reduced shear g. In the weak-lensing 
regime, characterised by <K 1 and 71,2 1, the ellipticity ap- 
proximates the shear, e ~ y. In gravitational lenses capable of 
strong lensing, the weak-lensing approximation typically fails 
very close to the centre. While the linear relation between ellip- 
ticity and shear can then be applied in the outskirts of the lens, 
corrections may become necessary near the core. 

Critical curves 9^ are closed curves in the lens plane con- 
sisting of points where cannot be inverted, det^(Oc) - 0. 
Their images in the source plane under the lens mapping (|5jl are 
called caustics; they are given by fic - 0^- Vi^(0c)- Obviously, 
the relation 



(8) 



is satisfied along critical curves. 

Any measurement of gravitational lensing which is based 
on local distortion information alone cannot distinguish be- 
tween lensing with the Jacobian matrices ^ and A^, with 
Ai^ 0. The matrix A^ produces images which are isotropically 
stretched by the factor A, but with identical ellipticity e as the 
matrix This invariance against the transformation ^ — > A^ 
causes k and y to be invariant against the transformation 



K ^ (I - A) + Ak , 7 — > /ly . 



(9) 



which obviously leaves the reduced shear Q invariant. This 
invariance tra nsformation was f i rst described as a "mass-s heet 
degeneracy" jpalco etalJll985t ISchneider & Seitjll995l) be- 
cause it tends to k ^ {1 - A) + Kin the limit of \1 - A\ ^ 0. 

We emphasise in the context of our joint reconstruction 
method that the critical curves are also invariant under the 
transformation Obviously, the defining condition det^ - 
for critical curves is unchanged if is multiplied hy A Q. 
Translated to the underlying lensing potential, the transforma- 
tion (|9j allows transformations of the potential of the form 

i/j ^ Aiff + ^—^{0^^ + ej) + a0i + bOi + c , (10) 



with arbitrary constants a, b and c. We shall later take advan- 
tage of this transformation for adjusting the lensing potential 
obtained from weak- and strong-lensing constraints. 

4. Outline of the method 

We propose to combine local constraints from weak and strong 
lensing on the lensing potential in the reconstruction of two- 
dimensional lensing mass distributions. We aim at the lensing 
potential because it is the smoothest lensing quantity available 
and because it provides a complete description of all lensing 
phenomena at least in the approximation of a single-lens plane. 
The lens is covered by a grid. The values i/r, of the lensing po- 
tential in the grid cells are considered as the model parameters 
and the method seeks to find a set of potential values {1/',) such 
that the agreement between the data and the model is optimised 
in a \t2&t-x^ sense. 

We introduce a function which is the sum of two terms, 
xl, and xl^ which encode information provided by weak and 
strong lensing, respectively. While x\, is defined on a grid cov- 
ering the cluster with low resolution, is restricted to a refined 



grid near the cluster core. Both contributions will be defined in 
the following subsections. 

4.1. Weak lensing 

Fol lowing the approa ch suggested bv lBartelmann et all (Il996l) 
and'S eitz et alJ ( Il998l) . one choice for the weak-lensing contri- 
bution to the x^ function is 



(11) 



where n is the number of grid cells covering the lens plane and 
e, is the expectation value for the ellipticity averaged within the 
/th cell. As already stated, e,- depends on second derivatives of 
the deflection potential tfr. 

The number of cells, n, must be chosen so that the grid 
cells are large enough for reasonably accurate measurements 
of the averaged ellipticities e,- , and yet small enough for the 
lensing potential iff not to change appreciably across a grid 
cell. We adopt the common assumption that intrinsic source 
ellipticities are randomly oriented and thus tend to zero when 
averaged within sufficiently large samples. The typical stan- 
dard deviation of int rinsic ellipticities from zero is cr^ ^ 0.3 
dBrainerd et al.ll996 ). Requiring that the noise of a shear mea- 
surement within a grid cell due to the intrinsic ellipticities be 
at or below the 10%-level, we need to locally average over a 
number of galaxies determined by 



0.3 



< 0.1 



(12) 



thus > 10. At a surfa ce density of, say, 30 arcmin" fo r the 
background sources (e.g. lFord et al .' 19%: R ix et al.l2004|) . this 
implies grid cells of > 35" side length. In the "concordance" 
ACDM cosmology, 1 /; 'Mpc spans ~ 5' at redshift ~ 0.3. The 
virial diameter of a massive cluster at that redshift is thus con- 
tained in fields of ~ (10' x 10'), approximately correspond- 
ing to ~ (20 X 20) pixels. This illustrates the grid resolution 
which we can expect weak-lensing cluster mass reconstructions 
to achieve. Then, the typical noise level c^^^,,■ per pixel is x 0.1, 
as follows from (I12> . 

The expectation value for ellipticities, e,-, depends on a com- 
bination of both the convergence and the shear. 



= I = 



where 1 - a: - |y| > (a) 



1-T^ elsewhere 



(13) 



(b) 



Expression (b) is restricted to the innermost region of the lens 
where both the convergence and the shear must be comparably 
large. 

Based on Eq. (I13> . the minimisation of xl, in dl li yields 
a non-linear relation between the measured ellipticities and 
the deflection potential. This technical proble m is conve- 
niently solved by an iterative procedure (see e.g.'Bradac e t alJ 
(.2005b.) "). Starting from /cf" = on all grid cells, subsequent 
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iterative approximations to the lensing potential are ob- 
tained by minimising 



(14) 



where g has to be chosen from two cases identified in ( I13> . 

Although this iteration provides an adequate solution for 
the potential, it may not be needed in actual reconstructions 
because of the difference in scales between weak and strong 
lensing. Strong lensing occurs typically within 30" to 1' from 
cluster cores, while the resolution that can be achieved by weak 
lensing is > 35". Thus, strong lensing is confined to the inner- 
most few cells of the weak-lensing grid. 

Once the potential values are found which reproduce the 
weak-lensing data, strong-lensing constraints can be added as 
described in the next subsection. 

So far, we have essentially reviewed the cluster recon- 
struction approach proposed by Bartelmann et al. (1994 and 
extend ed by maximum-entropy r egular! sation in i Seitz et al. 
( Il998h . While it was suggested in iBartelmann et alj C19 96) to 
minimise with the conjugate-gradient method, the finite dif- 
ferencing on the lensing-potential grid needed to obtain expec- 
tations for y, implies that, in the weak-lensing regime. 



= 



(15) 



is a linear equation in the i^, which can be solved using matrix 
inversion. The i^, are thus obtained from 



(16) 



with a sparse matrix S and a data vector 'V which are detailed 
in the Appendix. 

This was noted by iBradac et alJ ll2005albh . who have re- 
cently proposed and applied an alternative algorithm for com- 
bining weak and strong cluster lensing. The main difference 
between their and our algorithm is the way ho w constraints 
from str ong lensing are taken into account. While iBradac et all 
(|20p5b') use information from multiply-imaged sources, we 
advocate constraining the potential using the (approximate) 
knowledge of the critical curves from the location of large arcs, 
as described i n the next subsection. 

IPiego et'al. ( 2005 J recently proposed another cluster re- 
construction procedure which also combines weak and strong- 
lensing data. Their idea is to expand the projected cluster mass 
distribution into a set of basis functions which are then con- 
strained individually using shear and strong-lensing data. The 
minimisation then proceeds iteratively, using an adaptive grid 
to cover the cluster field. 



4.2. Strong lensing 

Strong lensing in clusters gives rise to highly distorted large 
arcs which occur in the immediate vicinity of critical curves in 
the lens plane. There, by definition (|8}, the positions of large 
arcs approximate the locations where det y[ = 0. If written in 
terms of the lensing potential i/', this condition translates into 
an expression which is quadratic in the second derivatives of (//. 



We suppose that a solution i/'^ has already been obtained 
from weak lensing using il6i . It exists on a coarse grid adapted 
to the low resolution of weak-lensing measurements. This grid 
is now refined near the cluster core to a resolution adapted 
to the critical curve(s). Since strong lensing typically occurs 
within the innermost arc minute around the cluster core, this 
grid will be much finer and smaller than the grid introduced 
for weak lensing. It thus forms a sub-grid which improves the 
resolution of the few central cells of the weak-lensing grid that 
contain the large arcs, and thus parts of the critical curve. 

In the idealised case of a fully known critical curve, we 
could identify a continuous chain of sub-grid cells covering the 
critical curve. Let the number of these cells be n*, then the con- 
tribution of strong lensing to the function could be 



X'. 
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(17) 



expressing the expectation that the Jacobian determinant be 
zero within the tolerance expressed by osi in all sub-grid pixels 
covering the critical curve. 

The derivative of the above function with respect to the po- 
tential values can be incorporated into the matrix approach as 
detailed in the Appendix. This leads a term to be introduced to 
the previous system il6\ . 



3j,if,,^^j-Tj 



(18) 



The implementation of this approach is addressed in the next 
subsection. 

The resolution of the central sub-grid introduced for the 
strong-lensing constraints has to be high enough to follow the 
critical curve with sufficient accuracy. The tolerance o-,, quan- 
tifies tolerable deviations of (det^l), from zero. It combines the 
uncertainty in the position of the critical curve, which must be 
estimated from the observed arc positions, and the deviation of 
det from zero expected within one sub-grid pixel. 

The second contribution can be suppressed to a negligible 
level because we are free to choose the sub-grid resolution. The 
first contribution can be estimated by considering the deviation 
of det^ from zero at distance 60 from the critical curve. 



OS 



de 



60: 



66 



(19) 



approximating the derivative of the Jacobian determinant at the 
critical curve by the inverse of the Einstein angle Oe- This is ex- 
act for the tangential arcs formed by an isothermal sphere, and 
reasonably accurate for similar lens models. Assuming uncer- 
tainties in the positions of critical curves of order 66 I" and 
Einstein radii of order 30", we find os ~ 3 x 10"^. 

We believe that our approach has two distinct advantages 
compared to methods using multiple images for constraining 
the lens model with strong-lensing data. First, it can be used 
simply based on the approximate knowledge of the location of 
the critical curves without the need to have or identify multi- 
ple arcs which are due to the same source. Prominent counter- 
arcs are often missing in strongly lensing clusters, which re- 
flects their lack of axial symmetry (.Grossman & Naravan.19881 
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ll989HKovne3ll989l) . Thus, methods relying on the identifica- 
tion of multiple images are only applicable to clusters which 
produce multiple large arcs, and in which the multiple images 
can be attributed to single sources with a high degree of cer- 
tainty. Second, methods based on multiple images require that 
all images identified as belonging to the same source be imaged 
on a single source. Even if an assumed model satisfies this re- 
quirement, it must be tested whether the source would have any 
additional images which are not observed. However, the inver- 
sion of lens models which is necessary to search for all images 
of a source is a highly nonlinear procedure. It thus appears to 
us that the combination of weak and strong-lensing information 
by identifying critical curves (or parts thereof) has substantial 
advantages compared to methods identifying multiple images. 

We emphasise that the method is fully local. Strong lensing 
constraints can be imposed for individual grid points. Although 
knowing the entire critical curves is desirable, only partial 
knowledge is necessary. Of course, the proposed reconstruction 
algorithm works best when the whole critical curve is available. 
Critical curves can be estimated from parametrised lens mod- 
el s. As an exa mple of the feasibility of this procedure, we refer 
to lSand etaP 1.2004,) and Meneghetti et al.. 12005.) . The authors 
show that observed arc positions constrain the whole critical 
line. 

In detail. iMeneghetti'et "ail (l2005l) perform a minimisa- 
tion constraining the critical curves that pass through the arcs. 
The reconstructed critical curves fit the real ones with high ac- 
curacy (see Sect. 5 for a quantitative comparison). It should be 
noted that such a method depends on the ellipticity of the de- 
flection potential. However, this aspect is not a weakness in our 
case because the weak-lensing reconstruction performed in the 
first step is able to supply that. 

Arcs in clusters may appear at different redshifts and thus 
trace different critical curves. Such information can be built 
into the;\f^ function by means of the redshift scaling t|2j. This 
could be achieved by identifying one critical curve (or arc sys- 
tem) as fiducial with a redshift Zs, and to scale the critical curves 
for all other systems according to their (spectroscopic or photo- 
metric) redshift by applying appropriate distance factors to the 
lensing potential i//. 

4.3. minimisation 

We now proceed to describe how the two contributions tox^ are 
joined. We begin by considering the weak-lensing constraint 
only, i.e. using the constraint leads to 

^ ^ ^ ^ = . (20) 
dtlfi 

This step may require an iteration if deviations from the weak- 
lensing limit need to be taken into account. 

The weak-lensing solution i/', obtained in this way on the 
coarse grid is improved on the fine grid as follows. Bi-cubic 
interpolation is used on the reconstructed potential near the 
cluster core in order to achieve a resolution high enough to 
incorporate the strong-lensing constraints. We emphasise that 
the interpolation is carried out on the lensing potential because 
it is the smoothest quantity available, and that the interpolation 



scheme has to be of sufficiently high order for the convergence 
and the shear to be continuous. In fact, since second derivatives 
of the potential are needed, interpolation schemes up to second 
order will fail. 

The location and shape of the fine grid onto which the po- 
tential is interpolated from the weak-lensing solution depends 
entirely on the strong-lensing constraints. There can be sin- 
gle or multiple arc systems which constrain the Jacobian de- 
terminant to zero at isolated locations near the cluster cores, 
or approximations to the entire critical curve may be avail- 
able from strong-lensing reconstructions of parametrised clus- 
ter mass models. 

We interpolate the potential onto a square-shaped grid en- 
closing all critical curves of our model cluster This is not only 
the simplest choice, but also driven by the idea that the central 
region needs to be described in closer detail. More complicated 
adaptive grids could be chosen to follow the critical curve(s) 
more specifically. 

On the fine grid, the interpolated weak-lensing solution iff* 
is improved by requiring 



= 



(21) 



which leads to the set of linear equations for ij/*, on the fine grid 



(22) 



where Ty is a data vector containing the information on the 
critical curves and the prime denotes that only few grid cells 
in the inner cluster region are constrained. Details are shown in 
the Appendix. If the interpolated weak-lensing solution already 
satisfies det^l - along the critical curve(s), T* - 0. Thus, 
T* quantifies by how much i/^^ needs to be corrected due to the 
strong-lensing constraints. 

We point out again that the potential reconstructed from 
weak and strong lensing alike allows transformations of the 
form 1I10> because neither the weakly-lensed ellipticities nor 
the critical curves are changed by them. This allows smoothly 
matching the potential values on the fine and coarse grids. 

5. Testing the method 

5.1. Synthetic data 

In order to test and illustrate the proposed method, we ap- 
ply it to synthetic background images lensed by a simulated 
galaxy cluster taken from the set described in Bartelmann et alJ 
1119981) . The cluster is simulated within the "concordance" 
ACDM model with a mass resolution of 1.0 x IO'^/i^'Mq. 
The current matter density parameter is ilmo - 0.3, the cos- 
mological constant is Qao - 0.7 and the Hubble parameter is 
h - 0.7. The cluster's redshift is z = 0.35, and its total mass is 
1.4 X IO'^/z^'Mq. It is embedded in a cube of 5/i^'Mpc (co- 
moving) side length, corresponding to 17.8 arc minutes. At a 
coarse-grid resolution of 32 x 32 pixels, one pixel has a side 
length of 33", in good agreement with the resolution constraint 
estimated for weak lensing above. 

Lensed data are synthesised by randomly distributing ellip- 
tical sources on the source plane, which we place at Zs = 1 for 
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simplicity. The sources have random intrinsic orientations, axis 
ratios which are drawn randomly from the interval [0.5, 1], and 
axes determined such that their area equals that of a circle with 
radius 0.5". The admitted simplicity of this choice of source 
properties should be irrelevant for demonstrating the feasibil- 
ity of our method. 

Starting from the known lensing potential of the cluster 
(projected along the three axes of the simulation cube), weakly 
and strongly lensed images are produced from the synthetic 
background sources by tracing rays through the deflection- 
angle field, which is the gradient of the potential. 

5.2. Results 

5.2.1 . Full knowledge of the critical curve 

Figure ^ shows the original and reconstructed radial conver- 
gence profiles. Far away from the cluster centre, the conver- 
gence is well recovered from weak-lensing data, but the central 
part suff'ers from the inevitable softening due to the resolution 
limit of weak lensing. This is corrected by adding the strong- 
lensing constraints. The agreement between the outer parts of 
the true profile and its weak-lensing reconstruction could be 
improved by means of the transformation ilO\ . but this would 
not remove the central discrepancy. 

The overall feasibility of the joint reconstruction is sup- 
ported by the agreement between the original critical lin es and 
those reconstructed using the algorithm by Meneghetti et alJ 
il2005h . Figure |2l shows the pixelised critical curve for one of 
the simulated clusters used to test the method, and its recon- 
struction. A direct comparison is possible because the simula- 
tion also provides the location of the entire critical curve. At 
the resolution used in our reconstruction algorithm, the agree- 
ment is completely satisfactory. The chain of pixels covering 
the critical curve is reproduced in great detail, providing ideal 
conditions for the constraint from minimising 

Figure shows two-dimensional reconstructions of two 
projections of the simulated cluster after applying the strong- 
lensing correction in their cores. The coarse-grid resolution far 
from the cluster centre is increased near the core where the 
critical curve of the cluster is located. The resolution of the fine 
grid is 16 times higher (per dimension) than that of the coarse 
grid, so that pixels of the fine grid have ~ 2" side length. 

The quality of the combined reconstruction is also illus- 
trated in Fig. 13 by a map of the relative differences between 
our reconstruction and the original convergence of the sim- 
ulated cluster. Evidently, the mass distribution is well repro- 
duced everywhere in the field (see the caption for a quantita- 
tive description), in agreement with the good recovery of the 
radial profile. Moreover, constraints from strong lensing sub- 
stantially improve the reconstruction in the cluster core. This 
demonstrates that the algorithm is working as expected. 

Figure|4]summarises the principal properties of our recon- 
struction method. The outer region (with low resolution) is well 
reproduced using only weak-lensing constraints, but the high- 
density peak near the cluster core is poorly resolved. The re- 
construction of the core structure is greatly improved by adding 
the constraints from strong lensing. Given the quality of the re- 



T ' ^ 

real profile 




comoving radius in Mpc/ln 



Fig. 1. Comparisons of the radial convergence profiles K(r) for 
two projections of the original cluster, their weak-lensing re- 
constructions, and the joint reconstructions after adding the 
strong-lensing constraints. While the profile obtained from 
weak lensing alone suffers from the inevitable softening, strong 
lensing considerably improves the agreement. 



suits, we did not need to apply the iterative algorithm suggested 
above for correcting the relation between measured ellipticities 
and shear. 

5.2.2. Approximate and partial knowledge of the 
critical curve 

Since our joint reconstruction algorithm rests on the approxi- 
mate knowledge of at least parts of the critical curve, we now 
investigate two main sources of uncertainty, namely position 
errors in the critical curve and incomplete knowledge of its lo- 
cation. For this purpose, we run reconstructions taking these 
aspects into account. Throughout, we do not assume precise 
knowledge of the simulated c luster's critical curve, b ut recon- 
struct it using the method by (^Me neghetti et alJl2005l see also 
ISand et al. 20 04). Positional errors of < 0.2" are found there. 
Since this is approximately the spatial resolution of the fine 
grid on which we reconstruct the convergence, we perform a 
reconstruction after randomly dislocating critical points by one 
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Fig. 2. Left panel: original critical curve of a simulated clus- 
ter. Right pa nel: critical curve reconst ructed with the algorithm 
proposed bv 'Meneghe tti et all ( l2005l) : see the text for details. 
Only the central part of the cluster is shown. The correspond- 
ing side length is a; 1' with pixel size of ^ 2". 




Fig. 3. Two simulated clusters are shown in the two rows. Their 
original convergence fields are displayed in the right column, 
their joint reconstructions using weak and strong lensing in the 
left. The colour scale is linear, the grey scales range from to 
1.1, and the contours are spaced by A/e = 0.1. 



pixel each. This procedure yields an upper limit to the effect of 
positional uncertainties in the critical curves. 

The left panel of Fig. |5] displays the relative difference 
between the convergences reconstructed from the randomly 
shifted and the original critical curves. The figure suggests that 
the algorithm is working well because the inaccuracies in the 
critical curve cause relative deviations at a typical level of a few 
per cent. 

Addressing the second source of uncertainty, we assume 
that only part of the critical curve is available (e.g. from the po- 
sitions and orientations of one or more arcs). The right panel of 
Fig- 13 shows the relative differences between one reconstruc- 
tion using the entire critical curve, and one using only 16% or 
15 pixels of it. Mimicking real observations, we further assume 
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Fig. 4. Deviation of the reconstructed from the original conver- 
gence fields, K and k, respectively. To avoid divergences near 
the field boundaries where both k and k are small, the grey 
scale encodes the relative difference between (1 -H/c) and (1 -i-^), 
i.e. (k-k)I(2+k+k). The grey scale ranges from -0.1 (white) to 
-hO.I (black). Obviously, the relative deviations are very small. 

that this fraction of the critical curve is split into two approxi- 
mately equal sections at the far ends of the critical curve. This 
approach simulates the presence of two arcs on opposite sides 
of a cluster. The relative differences are remarkably low, rang- 
ing around a few per cent and approaching 10% only in few 
points. 

The above tests confirm the underlying idea that a recon- 
struction based on joining weak and strong lensing features 
by using even part of the critical curve is feasible. The main 
sources of uncertainty do not substantially affect the method. 

6. Conclusions 

We have proposed a novel method for galaxy-cluster recon- 
struction which combines weak and strong-lensing data. The 
meth od is based on a lea st-y^ fitting of th e lensing potential if/ 
( Bartelmann et al. 1996.:.Seitz et alJl998l) and exploits the fact 
that the^^ minimisat ion can be carried oiit e fficiently by invert- 
ing a sparse matrix ( Bradac et al. 2005a b). Contrary to other 
method s proposed for joining weak and stron g lensing infor- 
mation dBradac et al.l2005btlDiego et all2005h . we propose to 
constrain the lensing potential obtained from the weak-lensing 
data by the approximate location of the critical curves, where 
the Jacobian determinant of the lens mapping must be close to 
zero. 

The angular resolution which can be achieved by weak- 
lensing cluster reconstructions is typically of order 0.5' for 
background source densities near (30 - 40) arcmin"^. This is 
much too coarse for tracing critical curves, whose typical sizes 



8 



M. Cacciato et al.: Combining weak and strong lensing in cluster potential reconstruction 



Fig. 5. Both panels show the relative deviations between two 
reconstructed convergences k in the centre of a galaxy cluster. 
As in Fig. 13 the grey scale ranging from -0.15 (black) to +0.1 
(white) encodes (k - k)/(2 + k + k) to avoid divergences for 
|a:| ^ 0. In both case s, the critical curve reco nstructed with the 
method proposed bv lMeneehetti et al1 ( l2005l) is used instead of 
the precise critical curve. Left panel: Critical points are ran- 
domly shifted in the radial direction by one pixel or ~ 2" prior 
to the reconstruction. While this procedure exaggerates realis- 
tic uncertainties, it provides an upper limit for the reconstruc- 
tion errors which are < 10% almost everywhere in the field. 
Right panel: Only 16% or 15 pixels of the critical curve are 
used, split into two sections opposing each other at the far ends 
of the critical curve. A random uncertainty in their positions is 
also included. Deviations are typically below ~ 10%. 



are (0.5 - 1)'. Thus, we propose to cover those regions of 
the cluster fields with a refined grid where strong-lensing con- 
straints are available. 

We test the performance of the method on synthetic images 
produced with simulated lensing clusters. We conclude that the 
mass distribution in galaxy clusters is well reproduced across 
the entire field, in particular where constraints from strong lens- 
ing features are introduced. In practice, our algorithm first finds 
a lensing potential on the coarse grid which fits the weak- 
lensing data best. This approximation to the potential is then 
interpolated into the cells of the fine grid covering the critical 
curves or parts thereof, and refined by ^lx^ minimisation taking 
the strong-lensing constraints into account. 

We have used several approximations here. First, we treat 
weak lensing to first order in the shear y, i.e. we compare 
measured ellipticities to y rather than the reduced shear g - 
y(l - /<•) '. Retaining the desirable linearity of the method, this 
can easily be overcome by introducing an iteration scheme in 
which the convergence k reconstructed in the previous step is 
used to update the r educed shear of the current step (see also 
iBradac et al .l2005hl) . Given the quality of our simulated results, 
we did not need to use this iteration scheme. 

Second, we have assumed knowledge of the entire criti- 
cal curve, which is of course not directly observable. Cluster 
images such as those recently taken with ACS on-board 
HST, however, show so many large arcs that critical curves 
can in fact be well constrained all around the cluster cores. 
Furthermore, critical curves can reliably be infeiTed from 
parametrised strong-lensing models, in particular when com- 



bined with dynamic al const raints from central cluster galaxies 
(e.g. Meneghetti et alj2005l) . 

Third, we have assumed the strongly-lensed sources to be 
all at the same redshift. If multiple arc systems at different red- 
shifts are observed, a lensing potential can still be reconstructed 
for one fiducial source redshift Is by applying the distance cor- 
rection factors defined in 

Thus, we believe that the simplifications used here are not 
at all restrictive, and that the method suggested here provides a 
useful alternative to those proposed earlier. Our tests with ran- 
domly displaced critical points and partial knowledge of crit- 
ical curves demonstrate that the method works well even in 
presence of positional uncertainties and gaps in the known crit- 
ical curve. 
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Appendix A: Linear equations ^or minimisation 

We summarise here some technical aspects of how the physical quan- 
tities K and y are related to the discretised deflection potential i/' in our 
approach. 

Derivatives of ijj are replaced by common finite-differencing 
schemes. We use 9 grid points for k, 7 grid points for 72 and 4 points 
for 7i. This allows any lensing quantity to be written as the multi- 
plication of a well-defined matrix with the vector of lensing-potential 
values, 



K, = , y] = gjj,f>j and 7? = gf^^^ 



(A.l) 



The matrices and are very sparse becaus e the finite- 
differencing schemes use only near neighbours (cf. iBradac et alJ 
l2005bl) . Algorithms exist for eflicient inversion of such matrices. 

Based on the finite-difl"erencing schemes expressed by the matrix 
Eqs.. lA.lk the minimisation of is reduced to a linear algebraic 
equation. Starting from the;^'^ for weak lensing, we have 



" I 



11 ^ 



= Z ^ [(^' - + - jMM 

i=\ »'f 

i= 1 



which can clearly be written in the form 



(A.3) 
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with the matrix 



(A.4) 



and the data vector 



(A.5) 



Similarly evaluating the constraints from strong lensing yields 
(det 



/I* 

z 

1=1 



1=1 
4det:?I, 

a- 



(1 - K.) 



dKi dy\ dyj 
dij/'., dif/'j, dif/'j, 



-4detj7I, 



(A.6) 



where the and the y. ' are obtained from the interpolated weak- 
lensing solution. 

On the refined grid, the;^^^ minimisation <A.3> is modified by 



(A.7) 



Obviously, if the weak-lensing solution already satisfies det ^ = on 
the critical curves, T = there, and no correction is necessary. 
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